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A METHOD FOR THE DESIGN OF SWEPTBACK WINGS WARPED TO PRODUCE SPECIFIED 

FLIGHT CHARACTERISTICS AT SUPERSONIC SPEEDS ^ 

By Wabrbn a. Ttjckbb 


SUMMARY 

One of the problems connected with, the sweptback wing is the 
difficulty of controlling the location of the center of pressure and 
hence the pitching moment. A method is presented for design- 
ing a wing to be self-trimming at a given set of flight conditions. 
Concurrently, the spanwise distribution of load on the wing is 
made to be approximately elliptical, in an effiort to maintain low 
wing drag. 

These flight characteristics are achieved by warping the wing 
out of a plane. The required warp is determined by the values 
of the coefficients of a four-term series describing the pressure 
distribution; these values in turn are determined from four con- 
ditions on the lift, pitching moment, and spanwise load dis- 
tribution. 

The method is directly applicable to several wing plan forms, 
including the triangle and the sweptback plan form with finite 
tips, under the restriction that the leading edge must be subsonic 
and the trailing edge muet be supersonic. The application to 
any specific problem is simplified to a rontine computational 
procedure by the presentation of certain basic data in tabular 
form. A discussion is given of some points to be considered in 
the application of the method to a practical case, and several 
representative examples are worked out. The resulting wings 
are shown to be ones which might practicably be built. 

INTRODUCTION 

The evolution of the sweptback wing for efl&cient flight at 
supei'sonic speeds has reached the point where the stability 
and control problems are being investigated. This situation 
implies that not only the lift and drag of the wing but also 
the pitching moment must be considered in relation to the 
au’plane as a whole. 

In order to be truly efficient at the design Mach number, 
the wing should produce the design lift coefficient without 
creating about the airplane center of gravity a pitching 
moment that would require a large deflection of the trimming 
device (with a correspondingly large drag) . In addition, it 
is generally desirable that tne spanwise distribution of lift be 
as nearly elliptical as possible and that any adverse pressure 
gradients on the wing be small so as to retard separation of 
the flow. These two conditions are not sufficient to guaran- 
tee that the wing drag wfll be a minimum because at super- 
sonic speeds the drag due to lift is also dependent on the 
chordwise loading; they are, however, conducive to low wing 
di-ag. 

The use of wings warped to produce a constant pressure 
over the surface has been proposed to eliminate the large 
adverse pressure gradients encountered with the flat wing. 


For a given plan form, however, a uniform pressure distribu- 
tion allows no control over the pitching moment. The wing 
waiq) necessary to produce certain other pressure distribu- 
tions has been derived (ref. 1), but these distributions do not 
lend themselves readily to the control of pitching moment; 
in fact, the conical nature of the pressure distributions fixes 
the center of pressure at the center of area for triangular 
wings. 

In the present report, data are presented from which the 
wing warp necessary to produce a certain type of pressure 
distribution may be determined. A development is then 
given in which certain constants appearing in the expression 
for the pressure distribution' are determined by conditions 
on the lift, pitching moment, and spanwise load distribution. 
In this manner a method is derived for designing a wing of 
given plan form, operating at a given supersonic Mach num- 
ber, to have a specified lift coefficient, a specified center of 
pressm’e, and a nearly elliptical spanwise load distribution. 
Although the pressure gradients are not controlled directly 
in the method, the type of pressure distribution used insures 
that for most reasonable design conditions the gradients wfll 
not be excessive. There is no reason to beheve that a con- 
figuration using a self-tr immin g device designed by this 
method will necessarily have a lower drag than will a simi- 
lar configuration using a flat wing and a deflected trimming 
device. The possibility does exist, however, and should 
probably be investigated. 

The method is apphcable to a wide class of wing plan forms 
shown in figure 1; the principal requirement is that the lead- 
ing edge must be subsonic and the trailing edge must be 
supersonic. The presentation is made in a form suitable 
for engineering use, and a table and computational form are 
provided so that the apphcation of the method is reduced to 
routine computation. 



FiairaB 1. — Plan forms to which the method is applicable. 
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I Supersedes recently declasslBed NACA UM LS1F03, 195U 
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SYMBOLS 


<r(l-X) 

=(l-*)[l-cra-X)] 




C 

Or 


Od 

Cl 

a 

c„ 


k=— 


local chord 
root chord 


mean 


, . , . 2 r 2(1+X+X*) 

aerodynamic chord, J (rdy= 


drag coefficient, 


local lift coefficient, 


Liftonchordwise atrip dy in width 
gc 


lift coefficient, 

gs 

pitching-moment coefficient, positive when pitch- 
ing moment tends to move wing apex up. 
Pitching moment 
gSe 


m 


^ .. Tip chord 

taper ratio, ^^^^c h ^rd 

cotangent of sweepback angle of leading edge 
(see fig. 2) 

mi cotangent of sweepback angle of trailing edge 

(see fig. 2) 

ilf free-stream Mach number 

•as /3m (see fig. 3) 

Pl static pressure on lower surface of wing 

fxt static pressure on upper surface of wing 

P lifting-pressure coefficient, 

g free-stream dynamic pressure 

rs— (see fig. 3) 
mx 


s 


wing semispan (see fig. 2) 


S wing area 

x,y rectangular coordinates parallel and normal, re- 

spectively, to free stream, with origin at wing 
apex (see fig. 2) 

x' distance behind leading edge measured in free- 

stream direction 

Xo distance of moment axis behind wing apex (see 

_ %• 5) 

x' distance of moment axis behind leading edge of 

mean aerodynamic chord (see fig. 6) 
z distance perpendicular to ary-plane, positive up 

ANALYSIS 

GENEBAL 

A convenient method is derived in reference 2 for finding 
the wing shape corresponding to a given pressure distribu- 
tion. In the present report, the lifting-pressure distribution 
over the wing is taken to he of the form 


J 

p=Ch'+o,'x+CAy\+o,Y 


( 1 ) 



where the axis system is that shown in figure 2, and (7/, Ck', 
Cz , and C^ are as yet arbitrary constants. Other terms 
could have been included in the series but the terms shoivn 
gave acceptable results without requiring undue labor. 
For convenience, the coefficients of the series may bo replaced 
by others similar in nature such that the lifting-pressure 
distiibution can be expressed by the folloiving equation: 


I C2 X ■ C3 

C.r'U'j/ \-\Cr.Cr'^ 




( 2 ) 


For purposes of calculation, the wing is assiuned to have 
no thickness so that the shape derived is actually the moan 
surface of the wing. Within the assumptions of the linear- 
ized theory, an arbitrary thickness distribution, symmetrical 
above and below the mean surface, can then be added with 
no effect on the lift and pitching moment. 

Suitable int^ations of the pressure distiibution over 
the plan form may be performed to obtain equations for 
the lift coefficient, the pitching-moment coefficient, and the 
spanwise load distribution. One condition may then bo 
imposed on the lift, one on the pitching moment, and two 
on the spanwise load distribution. This procedure results 
in four linear equations in the foiu unknowns GiICl, OjICl, 
OsIOl, and CJOi,. The values for these constants may then 
be substituted into equation (2), and the shape of the wing 
(that is, the warp) necessary to produce this pressure 
distribution can be found by the method of reference 2. 

The foregoing materia] has described the method in general 
terms. In the following sections more detailed descriptions 
are given: First, of the procedure used to find the warp 
corresponding to each component of the pressm'e distribu- 
tion; second, of the method used to determine the constants 
for the case of plan forms having pointed tips; and last, 
of the corresponding procedure for plan forms having finite 
tip chords. Although the determination of the constants 
is in principle the same for both types of plan forms, certain 
simplifications occurring for the pointed-tip case make not 
only the actual numerical work of determining the constants 
but also the exposition of the procedure simpler for this case. 
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WAHP 


The warps or wing shapes necessary to produce the several 
components of the pressure distribution given in equation 
(2) are first found separately as functions of the four constants 
Ci/Cx, CilOi,, CijCi,, and Ci/Ci;- Later, after numerical 
values of the constants have been determined by the condi- 
tions imposed on the lift, pitching moment, and spanwise 
load distribution, the separate shapes are superimposed 
to form the final warped wing. 

The general idea in finding each wing shape is first to 

determine the slope of the wing surface associated 

with the pressure distribution under consideration and then 
to integrate this slope in the r-direction to get the s-ordinate 
at any point (the direction of z is taken mutually perpendicu- 
lar to X and y, positive upwards). Of the available methods 
for finding the wing slope corresponding to a given pressure 
distribution, that presented in reference 2 was chosen for 
the particular problem. The principal advantage of this 
method is that it eliminates the need for considering z in 
the integrations involved and so simplifies the integrations. 

The slope of the wing surface corresponding to each term 
of equation (2) is found by apphcation of equations (8) and 
(17) of reference 2. The wing shape as a displacement from 
the z=0 plane is then found by integrating the slope in the 
r-direction; thus. 




dx 


( 3 ) 


The following equations result for the wing shapes corre- 
sponding to the four components of the pressure: 

xRz,^ — X ^ ^2Vl— TiV— 2 cosh“^ 






m 


m 


J. 

nr\ 


+ 


Vl— «^(l+r) cosh“* 




+ 


|7i(l-hr)l 

(1 -r) cosh-» 


22= 




L{VI=W- 


2r^cosh“' 




'•a(l+r) 


+ 


+ 


VI 




23= 








(4b) 


4.V 


TlV- 


i TiV^ cosh“’ 

a+ry+2q-n^(r+i^ 

2Vl^" 


+ 


cosh“^ 


1+nh- 




+ 


(l-r)»-2(l-7i^(r-r^ l-nh T 

2Vi=^^ 71(1 -r)J 




12—1071= 


4ir (. 3(1—71*) 


TiVyi — TiV- 6r* cosh“^ ~l+ 
3 ti*(1— 71=*) nr\ 


(1-712)8/2 [_• 


r6-97i=-|-27i^ 


(7^+r®)- 


2—371* 


(r— 7^— 


cosh~ 


l+7i*r 


+ 


7i(l-f-,r)| 

pe- 9 ti=*-1-271^ 2—371*/ ^ 

^ (t^-7^ 


(1-712)3/2 L 


(4d) 


The significance of the quantities r and n is most clearly 
seen by reference to figure 3. Calculations have been made 
of the quantities Ri, R%, R^, and Ri, which are in a certain 
sense the conical parts of the wing shapes, and the results 
are presented in figure 4 and table I. The figure is intended 
to be merely illustrative; the results in the table should be 
used for actual calculations. A study of the figure provides 
a quahtative idea of the various wing shapes. One interest- 
ing fact to notice is that no infinities occur at the center line 
(r=0) for the cases in which the pressure is proportional to 
X and to y*. (Compare with the shapes derived in ref. 1.) 



Figueb 3. — Definitions of r and n: r=—=—: 

Vi rnx’ Vi 

EVALUATION OF CONSTANTS FOR POINTED-TIP WINGS (X=0) 

The pressure distribution given by equation (2) is inte- 
grated over the plan form to obtain an expression for the lift 
of the wing. If the limits shown in figure 5 are used, the 
foUevraig equation expresses the value of the lift: 



( 5 ) 
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If the indicated operations are carried out and the lift 
coefficient is formed, the following equation results; 


, Oy , 2 — k Cj . 1 CIs . 1 Cj 


( 6 ) 


The pitching moment about an axis a distance Xo behind 
the wing apex may also he found by the following equation: 

f* dy f (r-ro) dx (7) 

2 Jo J elm '-‘L 



FroTTRE 4. — Values of if,, i?i, ifj, and ff^. 




Figube 4 — Continued. 



Figure 4 — Concluded 


Figure 4. — Continued. 
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TABLE I.— VALUES OF Ri, Ri, Rt, Rt 


r 

n 

0.1 



0.3 

Ri 

Rt 

Rt 

Rt 

Ri 

Rt 

Rt 

Ri 

Ri 

Rt 

Rt 

Rt 

n 

— CO 

-a 0820 

a > 


— 00 

-0.0870 

CD 

0.0250 

— GO 

-a 0930 

GO 

a 0240 


-0.4663 

-.0817 

o.mo 


-04275 

-.0866 

a 1099 

.0244 

-a 4810 

-.4810 

a 1076 

.0233 


-.3561 

-.0890 

.0579 

i/TB 

-.3623 

-.0860 

.0559 

.0228 

-.3707 

-.0926 

.0535 

.0218 


-.2917 

-.0798 

.0274 

TxSVM 

-.2979 

-.0848 

.0254 

.0208 

-.3064 

-.0916 

.0230 

.0107 

HB 

-.2463 

-.0784 

.0068 


-.2623 

-.0834 

.0047 

.0183 

-.2608 

-.0900 

.0023 

.0173 


-.2109 

-.0767 

-.0084 


-.2171 

-.0818 

-.0104 

.0161 

-.2255 

-.0884 

-.0129 

.0146 


-. 1825 

-.0749 

-.0198 

^^B iTTB 

-.1884 

-.0799 

-.0219 

.0127 

-.1969 

-.0865 

-.0244 

.0U6 


-. 1581 

-.0728 

-.0289 

.0106 

-. 1643 

-.0778 

—.0309 

.0097 

-.1729 

-.0844 

-.0334 

.0088 

.16 

-.1372 

-.0704 

-.0350 


-. 1436 

-.0765 

-.0379 

.0066 

1619 

-.0820 

-.0404 

'.0056 

.18 

1101 

-.0680 

-.0413 

.0045 

-.1248 

-.0780 

-.0438 

.0035 

-.1337 

-.0796 

-.(M58 

.0024 

.20 

-.1030 

-.0654 

-.0454 

.0014 

-.1091 

-.0704 

-.0476 

.0004- 

-.1176 

-.0709 

-.0501 

-.0007 

.24 

-.0770 

-.0597 

-.0509 

-.0046 

-.0816 

-.0648 

-.0630 

-.0056 

-.0899 

-.ons 

-.0566 

-.0067 

.28 

-.0525 

-.0538 

-.0533 

-.0100 

-.0680 

-.0587 

-.0654 

-.0111 

-.0669 

-.0652 

-.0680 

-.0123 

.32 

-.0332 

-.0474 

-.0533 

-.0149 

-.0393 

-.0623 

-.0556 

-.0169 

-.0476 

-.0587 

-.0682 

-.0172 

.36 

-.0167 

-.0407 

-.0516 

-.0189 

-.0228 

-.0456 

-.0638 

-.0200 

-.0310 

-.0620 

-.0566 

-.0213 

.40 

-.0029 

-.0340 

-.6488 

-.0220 

-.0085 

-.0387 

-.0507 

-.0233 

-.0166 

-.(W50 

-.0636 

-.0243 

.44 

.0068 

-.0263 

-.0442 

-.0243 

.0039 

-.0317 

-.0464 

-.0254 

-.0043 

-.0379 

-.0493 

-.0268 

.48 

.0205 

-.0108 

-.0390 

-.0266 

.0146 

-.0246 

-.0413 

-.0266 

.0067 

-.0308 

-.0443 

-.0281 

.53 

.0208 

-.0128 

-.0331 

-.0266 

.0239 

-.0181 

-.0354 

-.0208 

.0160 

-.0236 

-.0383 

-.0284 

.56 

.0376 

-.0059 

-.0268 

-.0247 

.0320 

-.0106 

-.0263 

-.0260 

.0242 

-.0166 

-.0321 

-.0276 

.60 

.0446 

.0009 

-.0201 

-.0228 

.0338 

-.0038 

-.0224 

-.0241 

.0311 

-.0097 

-.0265 

-.0263 

.65 

.0516 

.0089 

-.0116 

-.0192 

.0158 

.0044 

-. 0139 

-.0206 

.0383 

-.0014 

-.0169 

-.0222 

.70 

.0568 

.0164 

-.0029 

-.0142 

.0511 

.0120 

-.0063 

-.0166 

.0438 

.0064 

-.0084 

-.0173 

.76 

.0605 

.0232 

.0053 

-.0082 

.0660 

.0188 

.0030 

-.0096 

.0470 

.0134 

-.0001 

-.0114 

.80 

.0625 

.0289 

.0128 

-.0016 

.0573 

.0247 

.0105 

-.0029 

.0603 

.0104 

.0074 

-.0047 

.83 

.0029 

.0308 

.0156 

.0013 

.0677 

.0267 

.0132 

-.0002 

.0508 

.0215 


-.0020 

.84 

.0630 

.0325 

.0180 

.0040 

.0578 

.0284 

.0167 

.0026 

.0510 

.0234 


.0007 

.86 

.0628 

.0340 

.0203 

.0067 

.0675 

.0308 

.0179 

.0052 

.0509 

.0249 

^^BnTrSB 

.0034 

.88 

.0622 

.0351 

.0222 

.0092 

.0571 

.03U 

.0200 

.0078 

.0505 

.0282 


.0063 

.00 

.0613 

.0360 

.0238 

.0116 

.0563 

.0319 

.0216 

.0101 

.0498 

.0272 

TtTrfl 

.0083 

.91 

.0606 

.0361 

.0246 

.0126 

,0657 

.0322 

.0222 

.0112 

.0493 

.0276 

^^BiitisB 

.0094 

.92 

.0590 

.0362 

.0250 

.0136 

.0550 

.0323 

.0228 

.0122 

.0488 

.0277 

^^B iulTH 

.0103 

.03 

.0501 

.0361 

.0254 

.0144 

.0542 

.0323 

.0232 

.0130 

.0479 

.0277 

.0203 

.0113 

.04 

.0581 

.0360 

.0256 

.0162 

.0532 

.0322 

.0234 

.0138 

.0470 

.0277 

.0206 

.0120 

.05 

.0570 

.0356 

.0256 

.0167 

.0621 

.0319 

.0235 

.0144 

.0460 

.0275 

.0208 

.0126 

.96 

.0556 

.0351 

.0256 

.0161 

.0607 

.0313 

.0233 

.0146 

.0448 

.0270 

.0207 

.0130 

.07 

.05U 

.0343 

.0251 

.0162 

.0494 

.0307 

.0230 

.0148 

.0434 

.0264 

.0203 

.0132 

.88 

.0524 

.0331 

.0241 

.0159 

.0476 

.0296 

.0223 

.0146 

.0419 

.0264 

.0197 

.0130 

.09 

.0500 

.0316 

.0231 

.0152 

.0454 

.0281 

.0211 

.0140 

.0396 

.0240 

.0186 

.0124 

1.00 

.0460 

.0289 

.0207 

.6133 

.0421 

.0254 

.0187 

.0121 

.0359 

.0214 

.0163 

.0106 
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n 

r 


0.4 /' 



0.5 



0.6 
















Ri 

Hi 

Rt 

El 

Ri 

Pi 

Pi 

Pi 

Pi 

P, 

Pi 

Pi 




■** 

7 










— 03 


GO 

0.0229 

— CO 

-0. 1098 

CD 

a 0210 

— CD 

-a 1189 

CD 

a 0210 


-0.4011 

1010 

a 1049 

.0223 

-0,6026 

— . 1096 

a 1021 

.0212 

-0. 6163 

-.1186 

a0994 

.0203 

.04 

-.3800 

1003 

.0607 

.0207 

-.3924 

-. 1087 

.0481 

.0197 

-.4050 

-.1178 

.0453 

.0187 

• .06 

-.3166 

-.0991 

.0203 

.0187 

-.3282 

-. 1076 

.0177 

.0176 

-.3408 

-.1107 

.0147 

.0166 


-.2709 

-.0977 

-.0004 

.0162 

-.2824 

-.10© 

-.0032 

.0162 

-.2950 

-.1163 

-.0081 

.0142 

.10 

-=72367 

-.0960 

-.0166 

.0135 

-.24H 

-. 1046 

-.0184 

.0124 

-.2597 

-.1136 

-.0213 

.0114 

.12 

-.2070 

-.0941 

-.0271 

.0105 

-.2184 

-.1026 

-.0300 

.0095 

I-.2310 

-.1U0 

-.0330 

.0084 

.14 

1828 

-.0920 

-.0361 

.0075 

1943 

-.1006 

-.0390 

.0064 

-.2067 

-. 1096 

-.0421 

.0054 

.16 

-.1620 

-.0897 

-.0431 

.0044 

1740 

-.0979 

-.0450 

.0032 

-.1859 

-.ion 

-.0492 

.0023 

.18 

-. 1438 

-.0872 

-.0487 

.0012 

-.1652 

-.0958 

-.0610 

.0001 

-.1076 

-. 1040 

-.0«8 

-.0010 


-.1270 

-.0846 

-.0529 

-.0019 

-.1389 

-.0929 

-.0560 

-.0031 

1514 

1019 

-.0591 

-.0942 

.24 

-.0999 

-.0788 

-.0685 

-.0080 

-.im 

-.0872 

-. 0616 

-.0093 

-.1235 

-.0961 

-.0650 

-.0194 

.28 

-.0769 

-.0727 

-.0610 

-.0136 

-.0881 

-.0809 

-.0643 

-.0148 

-.1003 

-.0898 

-.0077 

-.0161 

.32 

-.0566 

-.0002 

-.0613 

-.0186 

-.0686 

-.0743 

-.0646 

-.0199 

-.0807 

-.0831 

-.0680 

-.0213 

.36 

-.0403 

-.0593 

-.0697 

-.0227 

-.0618 

-.0674 

-.0632 

-.0242 

-.0638 

-.0761 

-.0069 

-.0256 


-.0203 

-.0523 

-.0667 

-.0260 

-.0373 

-.0603 

-.0603 

-.0270 

-.0491 

-.0683 

-.0641 

-.0292 

.44 

-.0133 

-.0461 

-.0520 

-.0284 

-.0260 

-.0520 

-.0560 

-^.0303 

-.0363 

-.0614 

-.0601 

-.0318 

■/48 

-.0029 

- 0370 

-.0470 

-.0293 

-.0136 

-.0456 

-.0513 

-.0315 

-^.0250 

-.0539 

-.0654 

-.0333 

.52 

.0007 

-.0306 

-.0419 

-.0301 

-.0038 

-.0382 

-.0457 

-.0319 

-.0161 

-.0463 

-.9498 

-.0339 

.66 

.0149 

-.0234 

-.0356 

-.0294 

.0010 

-.0309 

-.0396 

-.0313 

-.0065 

-.0383 

-.0437 

-.0334 


.0203 

-.0171 

-.0295 

-.0277 

.0119 

-.0237 

-.0330 

—.0297 

.0019 

-.0318 

-.0333 

-.0318 

.05 

.0204 

-.0079 

-.0205 

-.01M2 

.0196 

-.0160 

-.0246 

-.0264 

.0090 

-.0224 

-.0289 

-.0283 


.0362 

.0000 

-.0120 

-.0194 

.0267 

-.0068 

-.0160 

-.0217 

.0154 

-.0139 

-.0204 

-.0242 

.75 

.0396 

^.0073 

-.0038 

-.0135 

.0303 

.0007 

-.0078 

-.0159 

.0204 

-.0060 

-.0121 

-.0186 


-^0422 

.0130 

.0038 

-.0069 

.0333 

.0074 

-.0001 

—.0096 

.0240 

.0010 

-.0044 

-.0121 

.82 

.0428 

.0158 

.0086 

-.0042 

.0341 

.0097 

.0027 

-.0067 

.0249 

.0036 

-.0014 

-.0095 

.84 

.0432 

.0177 

.0092 

-.0016 

.0347 

.0119 

..0063 

-.0040 

.0268 

.0053 

.0010 

-.0065 

.86 

-.0432 

.0195 

.0115 

.0012 

.0349 

.0137 

.0077 

-. 0012 

.0262 

.0079 

.0037 

-.0041 

.88 

.0431 

.0208 

.0136 

.0039 

,0349 

.0163 

.0098 

.0013 

.0265 

.0097 

.0058 

-.0013 


.0422 

.0218 

.0154 

.0060 

.0346 

.0165 

.0116 

.0037 

.0263 

.0112 

.0077 

.0010 

.91 

.0420 

.0223 

.0160 

.0073 

.0342 

.0170 

.0124 

.0050 

.0262 

.0117 

.0085 

.0024 

.02 

.0414 

.0226 

.0166 

.0083 

.0338 

.0174 

.0131 

.0069 

.0259 

.0123 

.0093 

.0034 

.03 

.0408 

.0228 

.0171 

.0092 

.0333 

1 rrB 

.0137 

.0070 

.0255 

.0127 

.0099 

.0044 

.94 

.0400 

.0228 

-0176 

.0100 

.,0320 

^^B 1 t??B 

.0141 

.0076 

.0251 

.0131 

.0105 

.^0063 

.05 

.0391 

.0227 

.0170 

.0106 

.0318 

^^B 1 T7!B 

.0142 

.0085 

.0244 

.0131 

.0103 

.0059 

.96 

.0380 

.0224 

.0176 

.OUO 

.0308 

^^B 1 vnfl 

.0143 

.0089 

.0237 

.0131 

.0109 

.0066 

,07 

.0367 

.0218 

.0173 

.0113 

.0298 

^^B 1 

.0141 

.0093 

..0228 

.0128 

.0103 

,0071 

.03 

.0351 

.0210 

.0168 

.0112 

.0283 

^^B 1 lisB 

.0137 

.0091 

.0216 

,0123 

.0105 

.0072 

.09 

.0332 

.0197 

.0158 

.0106 

.0266 

^^B~rT?fl 

.0128 

.0088 

.0200 

.0114 


.0068 

BB 

.0302 

.0173 

.0136 

.0090 

.0238 

Hill 

.0109 

.0071 

S.0176 

.0096 

.0081 

.0054 
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TABLE I— VALUES OF Ri, Ri, Rt, fii— Concluded 



n 

r 

0.7 

os 

o.g 


Si 

Ri 

Ri 

if . 

ifl 

Ri 

ifl 

ifl 

ifl 

ifl 

ifl 

if ) 

0 

— CO 

-0.1284 

CO 

0.0200 

OD 

-0.1384 

CD 

a 0192 

— 00 

-0.1487 

CO 

a 0184 

.02 ' 

-a 5288 

—.1282 

0,0068 

.0194 

^0.6431 

-.1381 

00941 

.0185 

-0.6582 

-.1483 

0.0911 

.0177 

.01 

—.4185 

-.1274 

.0425 

.0178 

-.4329 

—.1378 

.0397 

.0170 

-.4479 

-. 1476 

.0360 

.0162 

.06 

-.3541 

-.1262 

.0110 

.0167 

-.3686 

-.1882 

.0090 

.0149 

-.3836 

-.1464 

.0062 

0140 

.08 

-.3085 

-. 1248 

-.0090 

.0132 

-.3228 

—.1347 

-.0118 

.0123 

-.3378 

-. 1449 

-.0147 

.0115 

.10 

-.2732 

-. 1231 

-.0242 

.0104 

-.2876 

-.1330 

—.0272 

.0095 

.3000 

-. 1432 

-.0316 

.0087 

.12 

-.2446 

-. 1211 

-.0369 

.0075 

—.2687 

-^310 

-.0389 

.0065 

-.2737 

-. 1411 

-.0419 

.0060 

.14 • 

-.2202 

-. UOO 

-.0451 

.0043 

-.2346 

-.1289 

-. M 81 

.0034 

-.2494 

-.1389 

-.0551 

.0036 

.16 

-. 1904 

— . 1165 

-.0522 

.OOU 

—.2135 

—.1262 

-.0651 

.0001 

-.2284 

-. 1386 

-.0686 

—.0008 

.18 

-. 1812 

-.1140 

-.0579 

-.0021 

—.1962 

-.1238 

— .06 U 

-.0031 

-.2100 

-.1339 

-.0043 

-.0011 

.20 

-. 1647 

-.1113 

-.0623 

-.0053 

—.1788 

—.1210 

-.0666 

-.0064 

-. 1936 

-.1310 

-.0689 

-.0074 

.24 

-. 1368 

— .lOM 

-.0683 

-.0 U 6 

-.1608 

-.1160 

-.0717 

-.0128 

-.1654 

-.1260 

-.0752 

-.0139 

.28 

-. 1136 

-.0990 

—.0712 

-.0174 

—.1274 

-.1088 

-.0748 

-.0187 

-. 1419 

-.1184 

-.0784 

-.0200 

.32 

-.0937 

-.0922 

-.0718 

—.0227 

-.1076 

—.1017 

-.0756 

-.0241 

-.1218 

-.1 U 4 

-.0794 

-.0256 

.36 

—.0767 

—.0861 

-.0706 

-.0272 

—.9003 

-.0944 

-.0746 

-.0287 

-.1044 

-. 1040 

-.0786 

-.0304 

.40 

-.0618 

-.0777 

-.0680 

-.0309 

—.0752 

-.0869 

-.0721 

-.0326 

-.0891 

-.0963 

-.0764 

-.0343 

.44 

-.0488 

-.0701 

-.0642 

-.0336 

-.0619 

—.0792 

-.0686 

-.0353 

-.0766 

-.0884 

-.0729 

—.0373 

.48 

-.0374 

-.0625 

-.0506 

-.0353 

—.0500 

-.0713 

-.0637 

-.0375 

-.0637 

-.0804 

-.0686 

-.0394 

.62 

-.0273 

-.0547 

-.0541 

-.0360 


—.0634 

-.0587 

-.0382 

-.0530 

-.0722 

-.0034 

-.0405 

.56 

-.0183 

-.0470 

-.0482 

-.0357 

-.0307 

—.0666 

-.0529 

-.0380 

-.0434 

-.0641 

-.0677 

-.0405 

.60 

-.0106 

-.0304 

-.0418 

-.0343 

-.0225 

—.0478 

-.0468 

-.0389 

-.0350 

-.0660 

-.0616 

-.0396 

.65 

—.0022 

-.0301 

-.0336 

-.0313 

—.0137 

-.0380 

-.0383 

-.0341 

-.0261 

-.0458 

-.0433 

-.0373 

.70 

.0047 

-.0213 

-.0260 

-.0200 

-.0063 

-.0287 

-.0299 

-.0298 

-.0175 

-.0363 

-.0361 

-.0329 

.76 

.0102 

-.0130 

-.0167 

-.0214 

-.0002 

—.0200 

-.0216 

-.0244 

-.0108 

-.0270 

-.0266 

-.0275 

.80 

.0144 

-.0055 

-.0000 

-.0160 

.0046 

—.0120 

—.0137 

-.0181 

-.0062 

-.0186 

-.0187 

-.0213 

.82 

.0156 

-.0028 

-.0060 

-*'.0127 

.0062 

—.0091 

-.0108 

-.0154 

-.0033 

-.0163 

-.0166 

-.0180 

.84 

.0167 

-.0004 

-.0034 

-.0094 

.0076 

-.0063 

-.0079 

-.0124 

-.0016 

-.0122 

-.0126 

-.0158 

.86 

.0174 

.0020 

-.0007 

-.0068 

.0086 

-.0037 

-.0052 

-.0099 

-.0001 

-.0094 

-.0101 

-.0131 

.88 

.0180 

.0040 

.0016 

-.0041 

.0095 

—.0016 

-.0028 

-.0072 

.0013 

-.0068 

-.0072 

-.0101 

.00 

.0182 

.0058 

.0036 

-.0016 

.0101 

.0007 

-.0004 

-.0( M 6 

.0022 

-.0044 

-.0047 

-.0076 

.91 

.0181 

.0066 

.0046 

-.0006 

.0103 

.0015 

.0005 

-.0033 

.0027 

-.0032 

-.0035 

-.0062 

.02 

.0180 

.0072 

.0054 

-0007 

.01 (H 

.0023 

.0014 

-.0021 

.0031 

-.0022 

-.0024 

-.0049 

.03 

.0170 

.0078 

.0061 

.0017 

.0104 

.0031 

.0022 

-.OOU 

.0034 

-.0013 

-.0017 

-.0033 

.04 

.0178 

.0082 

.0067 

.0028 

.0104 

.0037 

.0029 

.0000 

.0037 

-.0004 

-.0007 

-.0025 

.05 

.0171 

.0086 

.0072 

.0035 

.0102 

.0042 

.0035 

.0009 

.0038 

.0002 

-. 0001 

-.0017 

.06 

.0166 

.0087 

.0076 

.0041 

.0099 

.0046 

.0040 

.0017 

.0039 

.0009 

.0007 

-.0000 

.07 

.0158 

.0086 

.0076 

.0046 

.0095 

.OMS 

.0044 

.0023 

.0038 

.0013 

.0011 

.0001 

.08 

.0150 

.0083 

-0074 

.0040 

.0089 

.0048 

.0043 

.0028 

.0038 

.0017 

.0016 

.0007 

.09 

.0137 

.0077 

.0069 

.0049 

.0080 

.0< M 5 

.0042 

.0027 

.0033 

.0017 

.0016 

.0011 

LOO 

« 

.0 U 7 

.0002 

.0055 

.0037 

-0066 

.0038 

.0031 

.0021 

.0024 

.0010 

.0009 

.0007 



After the pitching-moment coefficient is formed, the following 
equation is obtained;. 


II 

pak) 
L2 Cr' 

2-jfc - 
2(l~k)_ 

+ 

1 

“1 ak) 
_2<V 

S-k "I 

8(1-*)J 



2—k xo 3— O: 


!i+Fn 

-*) _ 




Xo 


4-k -\Ct 
'20(l-it)J ^ 


( 8 ) 


The spanwise load distribution is found from the following 
integration: , 

CCi 


CtCl 


Cr J Cj; 


dx 


where c is the local chord and c» is the local lift coefficient. 
The integration has already been made in finding the lift. 
The following equation results: 


cci f Cl ^ 1 — k Cj\ (0\ 7_ t7j (?3'\ 

CrC^TKOj 2 GJ OJ'" 

fl-^kC,,0, 0,\, G,. 


( 0 ) 


For purposes of reference, the spanwise load distribution if 
elliptical would be given by. the following equation: 


CtGz, TT 


( 10 ) 


where X=0 for the pointed-tip case. Equations (6), (8), and 
(9) are now used to find values for GijOi,, GjIGl, Gz/Gl, and 
GijGi,. The following conditions are first applied to equation 
(9): 


' cc,\ 2(l-fX) 

PtGjj) ffmmO ■7T 


(value for ellipse) 





(value for ellipse) 


( 11 ) 


where again X=0 for the pointed-tip case. These conditions 
are not quite arbitrary hut were chosen after trial of a 
number of possibilities. The selection of these particular 
conditions not only made possible a solution for the four 
imknowns but also resulted in a single equation for the 
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spanwise load distribution which was a fair approximation 
to an eUipse. (The degree of the approximation is shown 
subsequently.) 

The solution for the constants may now be made. Sub- 
stitution of equations (11) into equation (9) gives the follow- 
ing values for Ci/Ct and CiIOl' 

Cl 4k . 1— i C 3 

Gj~{l+kWl+k Gi. 

Ci_ 4 ^ C 3 

Cz, l-j-A Ci, 


These values are substituted into equation ( 6 ) and the foUow- 
ing solution for GtJCi, is got: 


Gz 



When equations ( 12 ) are substituted into equation (9), the 
following equation for the spanwise load distribution is 
obtained; 

This equation is xmiformly valid for ah. pointed-tip wings, 
independent of the lift, center of pressure, Mach niunber, 
and relative sweeps of the leading and trailing edges. The 
load distribution given by equation (15) is compared with the 
elliptical distribution of equation (10) in figure 6 . As a 
matter of incidental interest, the spanwise center of load on 
one wing panel is located 0.409 semispan outboard of the wing 
center line for the load represented by equation (15) as 
compared with a corresponding value of 0.424 for the 
elliptical load. 


The values of CiIGl, CifCz, and CJCl are substituted into 
equation ( 8 ) and the left-hand side is set equal to zero to 
arrive at the solution for GsICl. The solutions are collected 
in the following equations: 



03 

8(l+i) 

— 1 
CO 

1 

r 

1+7*+*“ 

3 Xq 

Gr.' 

l—k 

Lio(i-*) 

5a-(l— *“) 

23 ;. 

Qz 

8(1+*) [ 

" 7 - 3 * 

1+7*+*“ 

Xi~\ 

Cl 

" 1-* 1 

-10(1-*) 

5ir(l-*“) 

C J 

c. 

4 

2 Gi 



Cl 

(l+*)7r 

l+kCL 



Cl 

4* 

, 1-* <73 



Cl 

(l+*)ir 

^1+*^ 




The two forms for C^ICl are given because in some cases the 
center of pressure is located more conveniently with respect 
to the mean aerodynamic chord c. The geometrical relation 
between Xo and W is indicated in figure 5; the analytical 
relation is 


Xq_ 2 X' . 1 

<v“3 c ' 3(1—*) 


(13) 


For k—Q (triangular wing), equations ( 12 ) simplify to the 
following equations : 


or 



^3^0/6 1 3 rp\ 

Vz \5 5ir 2 cj 

<^ 3 ^ 0/7 1 ^ 

Cjr, \10 5 tt cJ 

Qi 4 „ Gi 


(14) 



Figubb 6. — Spanwise load distribution for pointed-tip wings compared 
with elliptical load distribution. 


BVALttATION OF CONSTANTS FOE WINGS WITH FINITE TAPEE 

The problem of wings with finite taper can be approached 
in two ways. The more obvious method is to assume that 
the pressure distribution defined by equation ( 2 ) appHes 
over the entire wing surface and to calculate the required 
warp. (A separate calculation for the warp in the tip region, 
which is shown shaded in fig. 2 , is required, but this calcula- 
tion is not impossible to make.) The disadvantage of this 
procedure is that at the tip the required wing slope takes on 
very lai^e values (theoretically infinite). A more practical 
approach, and the one adopted in this report, is to relax the 
condition on the pressure in the tip region. For a flat lifting 
wing with subsonic leading edges, the average pressure in 
the tip region is known to be dose to zero (refs. 3 and 4). 
It is not entirely illogical to suppose that for a shghtly warped 
wing the pressure in the tip will also be very small. If for 
the warped wing the pressure in the tip were taken to be 
exactly zero, the equations for the lift, pitching moment, and 
so forth would be derived by first integrating the pressure 
distribution defined by equation ( 2 ) over the entire wing, 
including the tip region, after which the integral of the same 
pressure over the tip region would be subtracted. In order 
to keep the equations within reasonable limits, a constant 
pressure was instead subtracted from the tip region. The 
value of this pressure was taken to be the value given by 
equation ( 2 ) at the middle of the tip chord; namdy. 


Ci_0, 

CzTCz 


Gz 


X(l- 

-*)■ 


1 ^3 

1 

‘ 2(1- 

-x)_ 





( 16 ) 
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After the foregoing assumption has been made, the de- 
termination of the four constants proceeds very similarly to 
that for the pointed-tip wing: The pressure distribution is 
integrated to obtain equations for the lift, pitching nioment, 
and spanwise load distribution; two conditions are set on the 
spanwise load distribution, and the resulting, two equations 
are solved together with the equations for the lift and the 
pitching moment to give values for the constants OiIOl, 
QiIOl, OzIOl, and CJCl. The procedure just indicated is 
now given. 

The hft equation, corresponding to equation (6) for the 


pointed-tip wing, is 


1 = 



2-k 
3(1 -XI 


X*(H-/:) 2-X(l+Z:) ."ICa, 
3(1 -X*) 2(1 -X) AW:/ 


where 


“1-1-2X 

_3(1+X) 


l(73.ri+3X 

^J^+Lfid+x) 

G, 

Gl 

(17) 

\^(l-k) 


(18) 

^-(k+n){l-\^ 



In thifi equation, aU the teims containing A arise from 
integration of the pressure over the tip region, whereas the 
remainder of the terms represent the mtegration of equation 
(2) over the whole wing area. If X is set equal to zero, the 
equation reduces to the form given in equation (6). 


The pitching-moment equation, corresponding to equation (8) for the pointed-tip wing, is 


where 


(H-X+X*) 


&-={D 

p 

[ 


3(H-X) 


,"|ro r3(l-A:X) a+*)(l-X)'' 

_ Cr \_2{l—k) 2(1— Ar) 




~k\) {i+k)Ci~\) 


(l-k)+2il-\) 
2(1 -X) 


l-i , 3it , F(l-X) (\+k+k^(l-\f 


{l-x“^2"^ 1-k 4(\-k) 


X(l— i)-f2(l— X) 
2(1 -X) 




k(l-\) 3(l+Ar)(l-xy 


8(1 -i) 




/l-}-3X A\ Xq 

ri . 3m-x) 3(H-*)(i-x)* /-\\ 

V 4 

12 ‘ 8(l-k) 10(1-^) J/ 


3X^(1 -A:) 

2{t-l-n)(l-X) 

3X* , X='(27i-ht)(l-*) 

^^2{k+n) ' 2(k+nf{l-X) ^ 


(19) 


(20) 


The remarks following equation (17), with A replaced by B and O and with equation (6) replaced by equation (8), 
also apply to equation (19). 

Because of the existence of a tip region for the wing of fini te, taper, the spanwise load distribution must be de- 
scribed by two equations, one apphcable from the center hne out to the beginning of the tip region and the other 
apphcahle over the remainder of the span: 


For O^cr^l 


X(l-Ar) 
(t-fn)(l— X)’ 


CCi 

CtGl 


rw , i-k c,\ r., r.O: Of] r(n-i)d-x) (7. 

\^+2(i=x) wr^ wrwi 2 


<r^-(l-X)^o® 


(21a 


For ^0-^1, 

(A:+n)(l— X) 

cc, r(l+«)(l-X) W , (l-^)“-[2-X(l+A:)][X(H-7i)-(fc-l-n)] <7* \-\{l+n)-(k+n)A f (1- f n)(l-X) <7, 

1-ifc W:/ 2(l-i)(l-X) L 1-* A\CLCjj \ 1-* Gu 

(l-X)[2A<l-A:)-2(^-|-7i)-j-X(l-hA:)(A:-h7i.)] G^ (l-A:)-(l-X)(Z:-hn) G^ (l-\)(Jc+n) G^^\ _ 

2(l-t)(l-X) Gl 1-k G:/ 1-k Gj,f 
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Tho conditions given in equations (11) are now applied to equation (21a) to obtain the following relations: 

(7i 4k(l+\) 1—k Os 

GiT (1+*)^ ■*'(1-X)(1+*) ^ 

(7i _ 4(l— ^ Os 

Oz, (l+i)ir 1+i Ol 


( 22 ) 

(23) 


These values of Ci/Cx and OsjOz are substituted into equation (17) and equation (19). The center of pressure is fixed 
at ro by putting 0„ equal to zero in equation (19) and the following equations for OsIOl and OJOl are obtained: 


O, 6(1 +X) r 8(1+X)-4X^ (4-2X)(l+X) ' 

Ol (1+3X)-6(1+X)A L 35t r 


(241 


(1-A)(1+X+X*+X^ Os , r/l+3X % n , 3i(l-X) 3(1+*)(1-X)» O, 

8(l-X)(l+i) 4 ^)cr L2'^4(l-i) 10(1-*) 


{[3(l+*)-(l+*)(l-X)*-(2-X)(l+*)B] 


r(3+*-**)(l +X)-(3+5*+*^X^+(l+*+*^X» 


L" 


2(1-*) 


-(2-X)(l+*)(7 


']} 


(25) 


where A, B, and 0 are defined in equations (18) and (20). 

Equations (24) and (26) can be solved for C,/Ci and 
OsIOl, after which OijOt and OsJOl can be found from equa- 
tions (22) and (23). This calculation is best done with the 
numerical values for X, *, n, and ro/c, for the particular wing 
under consideration. The procedure is illustrated by an 
example in a subsequent section entitled “Numerical Ex- 
amples.” The relation between ro/cr and F/c corresponding 
to equation (13) for the pointed-tip wing is 

2(l-fX-hX») F 1+X-2X^ 
c,“ 3(l-fX) c'*'3(l-l-X)(l-*) 

. The spanwise load distribution corresponding to equation 
(16) for the pointed-tip wing can be obtained by substituting 
equations (22) and (23) into equations (21a) and (21b). 
The substitution into equation (21b) produces only added 
complexity, and the result is not given herein. The sub- 
stitution into equation (21a), however, gives the following 

simplified equations for 0 ^ <r ^ 1 — (1^^ ‘ 


C,(Jl L ’’’ Gr. I Ox 


(27) 


Unlike the spanwise load distribution for the pointed-tip 
wing, the load distribution for the finite-taper wing cannot 
be compared with an elliptical distribution for all wings but 
must be compared separately for each example investigated 
because of the form of equations (27) and (21b). The 
elliptical distribution is stiU given by equation (10). 


NUMERICAL COMPUTATIONS 
DEVELOPMBNT OF FORM FOE COMPUTATION 

After numerical values have been found for the constants 
OiIOl, OsIOl, OsIOl, and OijOL, these values are used with 
equations (4), or rather with the numbers in table I com- 
puted from these equations, to find the z-displacement 
corresponding to each component of the pressiue distribu- 
tion. The four displacements are then added to produce the 
final shape of the warped wing. In principle this process is 
straightforward so that in practice it may be reduced to 
routine computation. A form suitable for such computation 
is now developed. The particular form presented is one such 
that, at a given spanwise station a, the s-ordinate as a 
fraction of the local chord c is given as a function of x'/c, the 
fractional distance behind the leading edge of the local chord. 
As a typical example, the z-ordinate corresponding to the 
second term of equation (2) is considered. Erom equation 
(4b) the following relation is written: 

m Zg „ Os mcr { a: W CrV c 

OlC^Ol s 

The following geometrical relations are easily verified: 

^ .n Yi -\\i ^ I °^(1 X) 

--[l-o(l-X)]--l- 


^=l-v(l-X) 

mCr 1—* 
s ~1— X 
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SO that 

m 32 „ <7s (1— Z;)[l— 0<1— X)] (a/ , <Kl— >^) 

C7z> c ^ ' Ci. 1-X ( c ■^(l-i)[l-cr(l-X)] 

Now if ^ is defined as 

^ (l-Z:)[l-cr(l-X)] 

then 


to; plot, for the particular wing under consideration, x'jc 
against r (for the various values of a) from equation (30). 
Only)-values of x'lc between 0 and 1 are of interest. This 
"Ta;ct determines, for each v, the range of r to be used, and 
from the values of r used as arguments in table I those 
which give a satisfactory distribution of points along the 
chord are chosen. Kelations similar to equation (29) for 
the other s-components can be derived; these relations are 



<r (1 — Ar)[l — <7(1 — X)] 


m 3i 


A 1-X 


Cl c~ 

and it can be shown that 


m 33 




Ol c 


A , v(l— X) x' . 

r c ‘ (1— *)[1 -ct(1-X)] c ‘ ^ 

(28) 

m Zi_ 

so that 



Cl c 


m 32 Tf Oa a- /AV 
Cj, c^'DIaXt) 

(29) 

The amount of warp 


Thus, at a particular spanwise station <r on a given wing, 
32/c is a fimction only of r (since B as a fimction of r is known 
from table I). From equation (28), x'jc is also a function 
of r; thus, 

c r 

so that Zilc can be calculated for various values of x'jc. By 
experience, it has been found that the most satisfactoiy 
way of choosing values of r to use in the computations is 


0\ .^0 
Ol r 


(31) 


<7i and inversely proportional to m (for a given n) . 

For some purposes, the wing shape can be more conven- 
iently expressed in terms of s/Cr and xjCr rather than in 
terms of s/c and a//c. This conversion is easily made by 
use of the following equations: 


J=[I-a(l-X)] 

Ct V 


l-A: 


(32) 
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The foregoing equations are embodied in table II, which is 
self-explanatory. This form has been used to compute 
several examples, some of which are discussed in the next 
section. 

NUMEEICAL EXAMPLES 

Example I. — ^The first example chosen has the following 
characteristics: 

n=0.6 

i=0.6 

X=0 


CL=0 at 


^=0.25 

c 


This set of characteristics represents a sweptback wing, 
tapered to a point at the tips, with the center of pressure a 
little more than 60 percent of the mean aerodynamic chord 
ahead of the location for the corresponding flat wing. The 
lift coeflficient, sweepback angle of the leading edge, and 
Mach number are not specifled; the final amount of warp 
is directly proportional to the lift coefficient, and any 
combination of sweepback angle and Mach number that 
gives 71=0.6 may be chosen. 

The four constants are found from equations (12) to have 
the following values: 

^=4.4530 

^=-19.0819 

^=15.9021 

^=0.9071 

Several values of <r are selected and, for these values, plots 
of x'jc against r are made from equation (30). Shde-nile 
accuracy is sufficient for these calculations, and only a few 


x7c 



Fiqubb 7. — Plot of x'lc against r for example I., 



Fiqttbb 8. — Wing sliape for example I ■with ordinates expressed as 
fractions of the local chord c. 

points need be taken to define the curves. The resulting 
curves for example I are shown in figure 7. These curves 
are used to pick values of r from table I. The corresponding 
values oi i?i, R^, and Bi from table I are entered in table 
n, together with the other necessary data, and the indicated 
computations are carried out. 

The results of the computations are shown in figure 8, in 
which the ordinates are given as fractions of the local chord 
and the origin of the axes is at the leading edge of the local 
chord for each value of o-. Several features of the wing are 
evident from this form of presentation: namely, the reflex 
curvature of the airfoil sections near the center of the wing 
(the angle of attack is infinite at the center line), the disap- 
pearance of this reflex cmwatme at outboard sections, the 
relative twist between inboard and outboard sections, and 
the (variable) dihedral. A better picture of the actual •wing 
is obtained by plottiog the results as in flgure 9. In order 
to give more physical meaning to the picture, the results 
have been plotted for a lift coefficient of 0.2 and a leading- 
edge sweep of 60° (m=0.677). This last value thus corre- 
sponds to a Mach number of 1.44 since n is equal to 0.6. 
There are two points worth mentioning with regard to flgure 
9. The first is that, within the accuracy of the linearized 
theory used in this report, an arbitrary z(tr) may be added 
to the vertical ordinates without changing the aerodynamic 
characteristics of the 'wing. As pointed out in reference 1, 
this procedme is permissible so long as the resulting ■wing 
does not he far from the 2=0 plane (that is, modification of 
the wing shape by addition of a set of ordinates which de- 
pends only on <r (not on a:) may be practiced in moderation). 
The practical significance of this point is that the wing shape 
may be modified by this procedure to simplify the problem 
of locating spars. The other point is that for most configura- 
tions the inboard stations of the wing, which are those having 
the laa^est warp, are buried within a fuselage and, therefore, 
present no structural problems. (The effect of the fuselage 
on the aerodynamic characteristics is discussed subsequently.) 




I.o B B .4 .2 0 0 .2’ 4 £ a LO 12 L4 1.6 L8 2.0 2.2 ^4 2.6 

o- x/Cr 

Fiqube 9. — Wing shiipe for example I with ordinates expressed as fractions of the root chord Cr- Cx,=0.2; m=0.677; Af=1.44. 
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Tiiis particular example and presumably others not too 
extreme should therefore be quite practical to build. 

The pressure distribution for this example is shown in 
figure 10. Because of the far-forward specified location of 
the center of pressure, part of the wing carries negative lift. 
The spanwise distribution of load is that shown already in 
figure 6. 

Example n, — The second illustrative example has the 
foUowing'characteristics : 

71=0.8 

k=0 

X=0 

O„=0 at ==0.60 
c 


The values k=0 and X=0 characterize, a triangular wing. 
The center of pressure is at the same point as the center of 
pressure of the corresponding flat wing. The purpose of the 
present design is to show the kind ofwarp that might produce 
a wing with essentially the same center of pressure and span- 
wise load distribution as the flat triangle but without the 
steep pressure gradients that are known to promote leading- 
edge separation on the flat triangle, at least at low Reynolds 
numbers. A constant-pressure triangular wing, of course, 
has the same center-of-pressure location as a flat triangle and 
has no adverse pressure gradients, but the spanwise load 
distribution of such a wing is triangular rather than elliptical. 

The method of computation is much the same as that used 
in the previous example. The principal difference is that 
equations (14), rather than equations (12), may be used to 



413672—57 31 
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find the following values for the four constants: 

^==1.0907 

^==-0.9082 

^=1.0907 

^=0.9071 

The results of computations made with tables I and II are 
presented in figure 11. These plots clearly depict a wing 
the main part of which is almost fiat and which has a turned- 
down leading edge, a small twist from root to tip, and almost 




constant dihedral angle along the span. For easeof manufac- 
ture, this dihedral can be removed without affecting the 
aerodynamic characteristics very much. (See the discussion 
of this point under example I.) The pressure distribution is 
shown in figure 12; the chordwise gradients of pressure are 
not laj^e. The spanwise load distribution is again given by 
figure 6. 

Example m. — The third example chosen is a wing with 
finite taper, characterized by the following conditions: 

n=0.7 

k=0.6 

X=0.4 

(7»=0 at y=0.25 

These values are substituted into equations (18), (20), and 
(26) to obtain the following values: 

^5= 0.82857 

Cr 

A= 0.05861 
5=0.12308 
<7=0.20986 

Substitution of these values into equation (24) gives 

^'=- 0.1814 

Substitution of this value into equation (25) gives 
^=4.7438 

The remaining constants are obtained by substitution of this 
value of fls/Cj^into equations (22) and (23); thus, 

^=-5.2614 

^=2.6451 

From tbiR point the method of calculation is the same as that 
used in the two previous examples: suitable values of r are 
chosen, and the form of table II is followed to arrive at values 
for the wing ordinates. - The resulting wing shape is show 
in figure 13 together with the pressure distribution; the 
results have been plotted for a lift coefficient of 0.4 and a 
leading-edge sweep of about 59° (m=0.6), corresponding to a 
Mach number of 1.54. The center-of-pressure location 
shown in the figure for the fiat wing was found from reference 
5. The wing shape is not extreme, and the previous remarks 
concerning Ihe removal of the dihedral angle apply equally 
well to this case so that the wing can be built feasibly. The 
pressure plot shows the result of the assumption regarding 
the pressure in the tip region. 

The following spanwise load distributions are found from 
equations (27) and (21b): 

For 0.795, 

CCt 


0 


1.0 


0.891-0.502ff*-h0.109o* 
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For 0.796 


-^=2.82-2.44o—0.49<r‘+0.11<r» 

CjOl 



/i t I t I I 1 1 1 1 1 1 

0 .2 .4 .6 .8 LO L2 L4 L6 L8 2j0 

X/Cr 

FiamiB 13. — 'Wing shape and pressure distribution for example IIL 
Ct=0.4; in=0.6; 1^=1. 54. 


the opposite wing panel. This condition is expressed by the 
following inequality: 

(34) 

In addition, because of the approximate nature of the 
assmnption regar ding the pressure in the tip region, cases in 
which the tip rfe^n^vers a large part of the wing should be 
viewed witji^paiition. 

ComputHig time. — The exact time required to compute a 
giveBu^Kample depends on such factors as whether X or k or 
bp'm^are equal to zero and the number of points taken to 
m^me the wing surface. The following time estimates are 
'^ven 4s representative of those required by using a manually 
Dperated calculating machine. To calculate the constants, 
K to 1 hour is required, and to calculate eight spanwise sta- 
tions, with 14 points along the chord at each station, 8 to 
12 hoitrs. 

Body effect. — In the derivation of the present method, the 
wing has been considered as isolated; whereas, in practice it 
is usually mounted on a body, on which may also be mounted 
a tail. The available information, both theoretical and 
experimental, is not yet sufficient to allow an accmrate 
quantitative prediction of the effect of the body for the 
general case. (See refs. 6 and 7 for a discussion of the prob- 
lem.) I Some qualitative estimates can be naade, and by 
referenlce to whatever experimental data may be available 
for cojifigurations resembling the particular example imder 
consideration rough quantitative corrections can be applied 
for the effect of the body. If the wing is mounted on the 
body so that the chord line at the juncture is, paraUel to the 
body center line, then ^6 lift of the combinatibn when the 
wing is at its design position with respect to the free stream 
will probably be dose to the sum of the lifts of ,th^ isolated 
wing and the isolated body. If, however, the'tw6^iu® con- 
nected so that when the whig is at its design position pi^h 
respect to the free streani the' body^ ’at zeio an^e of atiaij^, 


These equations are plotted in figure 14, which also shows an 
elliptical load distribution for comparison. The load distri- 
bution for the example being discussed is a fair approximation 
to the ellipse so that no large drag increase relative to the 
fiat wing is to be expected as a result of the specified forward 
location of the center of pressme. As a matter of fact, the 
spanwise load distribution of the fiat wing is itself not 
elliptical, so that the drag of the warped wing might well be 
less than that of the flat wing. 

NOTES ON PRACTICAL APPLICATION 

Range of applicahility. — ^The method described in the pre- 
ceding sections is directly applicable to wing plan forms of 
the types shown in figure 1. The locations of the various 
Mach lines shown in the figure relative to the leading and 
trailing edges and relative to Jhe center line are significant. 
The loading edge must be subsonic and the trailing edge must 
be supersonic; these conditions are expressed by the following 
inequality: 

(33) 



For the case of plan forms with finite tips (X?^0), the Mach 
line from the leading edge of one tip must not cross over to 


FiatiBH 14. — Spanwise load distribution for example m compared 
with eUiptical load distribution. 
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then the lift of the combination -will probably be somewhat 
less than the design lift of the wing. In all cases, the center 
of pressmre wiU probably be somewhat more rearward than 
that calculated when only tJie isolated components are con- 
sidered. If these foregoing statements, which are obviously 
conjectural in nature, are accepted, then some allowance 
can be made for the effect of the body by adjusting the 
design conditions of the wing. The body-interference prob- 
lem is neither different nor more serious for the warped wings 
considered in this report than for conventional flat wings, 
and all the preceding remarks apply equally well (or poorly, 
as the reader may judge for himself) to both types of wings. 

Off-design operation. — ^In the course of a flight, the wing 
may be required to fly at the design Mach number at atti- 
tudes other than that for which it was designed. Within the 
limits of the linearized theory used in the analysis, the prin- 
ciple of superposition applies. The lift (and pitching 
moment) of the warped wing at an attitude different from 
the design condition is therefore simply the design lift and 
pitching moment plus (or minus) the lift and pitching 
moment of a flat wing of the same plan form at an angle of 
attack equal to the angular deviation of the warped wing 
from its design attitude. When the wing is required to 
operate at Mach numbers other than the design value, how- 
ever, no simple method is available for estimating the change 
in aerodynamic characteristics, and even to calcvdate the 
properties by the use of the linearized theory is a practicably 
impossible job. An experimental test is the only way to find 
the answer. 

Applicability to other problems. — ^Although the derivation 
of the complete method has been limited to wings of the 
types shown in figme 1, with approximately elliptical span 
loads, the basic results presented in equation (4) and table I 
are applicable to other wings as well. For example, a 
derivation similar to that presented in this report could be 
made for sweptback wings with cross-stream tips, such as 
that shown in figure 9 of reference 3. It is also conceivable 
that in some cases the shape of the spanwise-load-distribu- 
tion curve might be determined by some condition other 
than that of low induced drag. The information presented 
in equations (4) and table I coifld be applied to such cases. 

As an example of an application of the basic data of equa- 
tions (4) to a problem of a type different from that discussed 
in the section entitled ‘TSTmnerical Examples”, the design of 
a triangiflar wing with approximately elliptical loading in 
both the spanwise and the chordwise directions is discussed. 
For convenience, this wing is called example IV. (In ref. 8, 
Jones has shown that, for a lifting surface of narrow propor- 
tions lying near the center of the Mach cone, the minimum 
value of the drag due to lift is achieved when both the span- 
wise and ihe chordwise loadings are .elliptical.) ^ 

The chordwise load distribution is found from integration 
of equation-(2) to be given by the following equation: 


QiIOl, and OJCi. are those of equations (14) : namely, 

7 ^ Cl 
Cl 4 - <?3 


Substitution of these values into equation (35) gives the 
following equation for the chordwise load distribution: 

Local lift <?5 a: , /4 3 , /o /rV 

Total lift 2^Acr/'^\ 3ffAcr/ 

The chordwise load is now specified to be zero at the trailing 

edge This procedure gives the foUo\ving value for 

a/CL. 

Cl 37t 

and the chordwise load becomes 




The load distribution given by equation (36) is compared 
with an ellipse in figure 16. The span\vise load distribution 
is also repeated from figure 6 for the sake of easy comparison. 

The wing shape is readily calculated from tables I and II 
and is shown in figure 16 for Cx,=0.2, M=1.2, and ti=0.3. 




(b) Spanwiae lo^. 

Figtthb 16 . — Load distribution for example IV. 


Local lift Cl X ,fCi ,1 6',^ /x-V > 1 Ci, fx\ 
Total lift“(7r oj W W 


(36) 


If, as in the previous examples, the conditions of equations 
(11) arq applied to equation (3)^ then thp spanwise load dis- 
tribution is given by equation (15), and the values of Gil C l, 


The drag coefficient for these .conditions has been found bj 
graphical integration to be approximately (7 d= 0.0081 
The drag coefficient of a flat wing at the same conditions it 
0.0091 if full leading edge suction is assumed or 0.0161 if nc 
leading-edge suction is assumed (no leading-edge suction ha. 
been assumed for the warped wing) . 
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FiatmB 16, — Wing shape for example IV. Ci,=0.2; Af=1.2; n=0.3. 

CONCLUDING REMARKS 

A method has been presented for designing a sweptback 
wing to ha've certain specified fiight characteristics at super- 
sonic speeds. For example, a -wing of given plan form, 
operating at a given supersonic Mach munber, may be 
designed to have a specified lift coefficient, a specified center 


of pressure, and a nearly elliptical span'wise load distribu- 
tion. As an aid in the calculations required for any specific 
case, certain basic data and a computational form are pre- 
sented as tables. The procedure is illustrated by sevei’al 
examples. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, V.a., May 11, 1951. 
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